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Abstract 

In  this  work  we  explored  the  selectivity  of  single  nucleobases  towards  adsorption  on  chiral 
single-wall  carbon  nanotubes  (SWCNTs)  by  density  functional  theory  calculations. 

Specifically,  the  adsorption  of  molecular  models  of  guanine  (G),  adenine  (A),  thymine  (T), 
and  cytosine  (C),  as  well  as  of  AT  and  GC  Watson-Crick  (WC)  base  pairs  on  chiral  SWCNT 
C(6,  5),  C(9,  1)  and  C(8,  3)  model  structures,  was  analyzed  in  detail.  The  importance  of 
correcting  the  exchange-correlation  functional  for  London  dispersion  was  clearly 
demonstrated,  yet  limitations  in  modeling  such  interactions  by  considering  the  SWCNT  as  a 
molecular  model  may  mask  subtle  effects  in  a  molecular-macroscopic  material  system.  The 
trend  in  the  calculated  adsorption  energies  of  the  nucleobases  on  same  diameter  C(6,  5)  and 
C(9,  1)  SWCNT  surfaces,  i.e.  G  >  A  >  T  >  C,  was  consistent  with  related  computations  and 
experimental  work  on  graphitic  surfaces,  however  contradicting  experimental  data  on  the 
adsorption  of  single-strand  short  homo-oligonucleotides  on  SWCNTs  that  demonstrated  a 
trend  ofG>C>A>T  (Albertorio  et  al  2009  Nanotechnology  20  395101).  A  possible  role 
of  electrostatic  interactions  in  this  case  was  partially  captured  by  applying  the  effective 
fragment  potential  method,  emphasizing  that  the  interplay  of  the  various  contributions  in 
modeling  nonbonded  interactions  is  complicated  by  theoretical  limitations.  Finally,  because 
the  calculated  adsorption  energies  for  Watson-Crick  base  pairs  have  shown  little  effect  upon 
adsorption  of  the  base  pair  farther  from  the  surface,  the  results  on  SWCNT  sorting  by  salmon 
genomic  DNA  could  be  indicative  of  partial  unfolding  of  the  double  helix  upon  adsorption  on 
the  SWCNT  surface. 

[S]  Online  supplementary  data  available  from  stacks. iop.org/Nano/23/165703/mmedia 
(Some  figures  may  appear  in  colour  only  in  the  online  journal) 


1.  Introduction 

The  ability  to  enrich  single-wall  carbon  nanotube  (SWCNT) 
samples  by  a  specific  chirality  is  important  in  designing  mate¬ 
rials  for  improved  nanoelectronics,  while  diameter  selectivity 


could  be  useful,  for  example,  because  of  practical  consider¬ 
ations  [1],  SWCNT  sorting  with  ssDNA  has  been  recently 
demonstrated  based  on  earlier  work  [2],  in  which  d(GT)„  («  = 
10-45)  were  used  for  low-resolution  metal/semiconductor 
and  diameter  carbon  nanotube  separation.  Selectivity  of 
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Guanine  (G)  Adenine  (A)  Thymine  (T)  Cytosine  (C) 

Figure  1.  The  nucleobases  (NBs)  considered  as  mimics  for  ssDNA  adsorption  on  SWCNTs. 


the  SWCNT  to  particular  ssDNA  sequences,  dependent  on 
chirality  and  diameter,  was  evident.  However,  the  selectivity 
in  adsorption,  even  for  relatively  small  diameter  SWCNTs, 
has  not  been  fully  explained  yet.  Although  this  question 
was  addressed  by  molecular  dynamics  (MD)  simulations  [3], 
application  of  an  empirical  potential  is  limited  because  of  a 
lack  in  the  ability  to  model  SWCNT  o-n  re-hybridization 
with  diameter  or  subtle  changes  in  chirality,  among  other 
shortcomings.  For  example,  Johnson  et  al  claimed  [3], 
based  on  empirical  force-field  molecular  dynamics  (MD) 
simulations,  that  nucleotide  base  (NB)-SWCNT  binding  is 
not  affected  by  the  chirality  of  the  carbon  nanotubes,  yet 
this  is  contrary  to  experimental  work,  in  which  metallic 
SWCNTs  were  noted  to  be  less  eluted  [2].  Although 
recently  Tu  et  al  [4]  showed  a  targeted  elution  of  metallic 
SWCNTs  by  a  trial-and-error  ssDNA  mutation  approach,  a 
fundamental  understanding  of  intrinsic  interactions  is  still 
elusive.  For  example,  a  metallic  C(7,  7)  SWCNT  was  eluted 
with  a  modification  in  TTATTATTATTATT  to  C,  however 
considering  a  larger  diameter  (9.50  versus  7.72  A)  than  that 
of  the  C(8,  3)  SWCNT,  thus  rendering  the  two-dimensional 
wrapping  hypothesis  unclear.  Note  that  it  was  shown  by  MD 
simulations  that  solvent  and  entropic  effects  are  minimal 
in  this  case  [3].  Modeling  of  the  intrinsic  interactions  is 
important  also  in  understanding  the  enrichment  obtained  by 
the  use  of  salmon  genomic  DNA  (SaDNA)  for  SWCNT 
separation,  as  demonstrated  by  Kim  et  al  [5];  SaDNA  is  an 
abundant  byproduct  of  the  fishing  industry  and  could  provide 
cost  savings  in  large-scale  production. 

Quantum  mechanical  calculations  provide  higher  ac¬ 
curacy  than  empirical  force-fields  in  modeling  intrinsic 
interaction  energies,  but  have  been  limited  so  far  to  zigzag  and 
armchair  SWCNT  models,  performed  primarily  at  levels  of 
theory  that  may  not  take  into  account  nonbonded  interactions 
adequately  [6-9].  In  this  work,  the  objective  is  to  gain  an  un¬ 
derstanding  by  density  functional  theory  (DFT)  calculations 
of  the  selectivity  of  single  NBs  (guanine  (G),  adenine  (A), 
thymine  (T),  and  cytosine  (C),  see  figure  1),  and  of  AT  and 
GC  Watson-Crick  (WC)  NB  pair  model  systems,  towards 
adsorption  on  chiral  SWCNT  model  compounds.  Chiral 
SWCNTs  which  are  known  to  be  enriched  experimentally  [2], 
i.e.  C(6,  5)  (chiral  angle  27°)  and  C(9,  1)  (5.2°)  with  the 
same  diameter,  as  well  as  C(8,  3)  (15.3°)  having  a  similar 
diameter,  were  analyzed,  which  have  not  been  explored  so  far 


theoretically.  In  this  context,  we  discuss  aspects  of  modeling 
molecule-inorganic  surface  interactions,  and  the  limitations 
encountered  in  the  approach  taken.  Although  a  direct 
comparison  of  the  results  with  experimentally  observed  trends 
in  sorting  SWCNTs  by  single-strand  hetero-oligonucleotides 
or  double-strand  SaDNA  is  difficult  to  discern,  the  calculated 
NB  adsorption  energies  enabled  a  qualitative  explanation 
for  homo-oligonucleotides,  forming  the  basis  for  further 
exploration. 

2.  Computational  methods 

2.1.  Methods 

The  frequency-dependent  interactions  between  finite  systems 
can  be  expressed  in  terms  of  the  Hamaker  constants  C„, 
i.e.  A E(R)  =  —  where  R  is  their  separation.  The 

first  term  ( C(, )  for  two  molecules  A  and  B  is  given 
by  Cf  =  ^  f0°°  duaA(iu)aB(iu),  where  the  polarizability 
tensors  are  evaluated  at  the  imaginary  frequency  i u.  Wu 
and  Yang  [10]  proposed  that  the  exchange-correlation 
functional  can  be  corrected  non-self-consistently  by  the  term 
£d is  =  where  fd(R)  is  a  damping  function,  C(, 

is  based  on  the  polarizability  of  the  constituent  atoms,  and 
R 6  defines  the  asymptotic  form  of  the  interaction.  It  is 
assumed  that  quantum  mechanical  calculations  are  adequate 
in  describing  all  except  dispersion  interactions,  dependent 
on  the  level  of  theory  employed.  Consequently,  numerous 
schemes  to  include  long-range  dispersion  in  DFT  functionals 
for  modeling  larger  systems  have  been  suggested  [11,  12]. 
Notably,  Grimme’s  B97-D  has  proven  quite  successful  in 
modeling  intermolecular  interactions  for  molecular  systems, 
and  is  robust  because  of  extensive  validation  [13]. 

However,  accurate  prediction  is  further  complicated 
in  modeling  nonbonded  molecular-surface  interactions. 
For  molecular  adsorption  on  a  solid  surface  at  a 
small  distance  the  Hamaker  constant  is  given  by  C  — 
^  /0°°  dna(i»)(^“'j~[)  [14],  where  e(i u)  is  the  frequency- 
dependent  dielectric  function  at  the  imaginary  frequency. 
This  was  shown  to  be  important  in  modeling  interactions 
between  proteins  and  Si  and  Si02  surfaces  [15],  where 
the  TDDFT  calculated  surface  dielectric  function  and  the 
polarizability  of  the  adsorbate  were  combined.  Although 
predictions  of  the  dielectric  response  of  SWCNTs  from  first 
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Table  1.  Nearest  distances  (A)  of  heavy  atoms  in  NBs  from  the 
SWCNT  surface. 


C-SWCNT 

N-SWCNT 

O-SWCNT 

G 

C(6,  5) 

3.02 

3.17 

3.57 

C(9,  1) 

3.11 

3.06 

3.40 

C(8,  3) 

3.00 

3.20 

3.60 

A 

C(6,  5) 

3.12 

3.12 

N/A 

C(9,  1) 

3.10 

3.23 

N/A 

C(8,  3) 

3.12 

3.12 

N/A 

T 

C(6,  5) 

3.11 

3.15 

3.24 

C(9,  1) 

3.18 

3.22 

3.29 

C(8,  3) 

3.12 

3.27 

3.24 

C 

C(6,  5) 

3.10 

3.23 

3.38 

C(9,  1) 

3.12 

3.06 

3.51 

C(8,  3) 

3.06 

3.26 

3.26 

principles  have  been  reported  [16,  17],  theoretical  work  is  still 
lacking.  This  could  explain,  in  part,  the  weaker  adsorption  of 
metallic  SWCNTs  that  are  less  eluted  by  ssDNA,  as  observed 
experimentally  [2].  In  seeking  qualitative  trends  for  NB 
adsorption  on  chiral  SWCNTs,  in  this  work  we  assume  that 
the  solid’s  polarizability  can  be  approximated  by  those  of  the 
constituent  atoms.  The  functional  B97-D  was  applied  in  our 
work,  and,  for  comparison,  the  EFP2  method  [18]  which  has 
not  been  applied  as  yet  to  model  NB-SWCNT  interactions. 
This  method  provides  a  better  understanding  because  the 
interaction  energy  is  determined  as  a  sum  of  electrostatic, 
Pauli  exchange,  polarization,  dispersion,  and  charge  transfer 
terms  consistently  after  first  principles  [18].  The  combination 
of  results  from  the  two  methods  provided  insight  into  the 
interplay  of  electrostatic  and  London  dispersion  interactions 
for  NB  adsorption  on  chiral  SWCNTs. 

2.2.  Computational  details 


To  ensure  that  no  artifacts  in  the  starting  configuration  of  NB 
adsorption  on  the  SWCNT  occurred,  simulations  at  300  K 
for  1  ps  were  carried  out  with  the  self-consistent  charge 
density  functional  tight-binding  (SCC-DFTB)  scheme  using 
an  empirical  correction  to  the  dispersion  (DFTB-D,  mio-01 
parameters)  [19-21]  and  an  Anderson  thermostat.  Structures 
were  then  optimized  with  B97-D/6-31G*  [13,  22],  Note  that 
the  translational  unit  cells  for  the  chiral  SWCNTs  contain  388 
and  364  carbon  atoms  for  C(8,  3),  and  C(6,  5)  and  C(9,  1), 
respectively,  too  large  for  feasible  computations,  and  therefore 
about  half  of  the  unit  cell  was  used,  specifically  167  and 
161  atoms,  respectively.  Calculations  with  EFP2/6-31G  were 
performed  for  C(5,  0)  because  of  the  system  size,  where  the 
starting  structures  were  optimized  at  the  B97 -D/6-3 1G*  level. 
To  mimic  dsDNA  adsorption,  calculations  for  WC  models 
of  base  pairs  (AT  and  GC)  and  of  double  base  pairs  (ATAT 
and  ATGC)  were  carried  out.  All  calculations  were  performed 
with  Gaussian09  [23]  except  for  EFP2  with  GAMESS  [24]. 


(a) 


Figure  2.  Structures  of  systems  with  NBs  adsorbed  on  (a)  C(6,  5), 
(b)  C(9,  1 )  and  (c)  C(8,  3)  SWCNTs. 


3.  Results  and  discussion 

The  geometrical  parameters  of  the  optimized  SWCNT-NB 
molecular  systems  are  listed  in  table  1.  The  results  show 
stabilized  stacking  of  the  NBs  on  the  SWCNT  surface 
(see  figure  2).  jr-stacking  of  DNA  NBs  on  surfaces  is 
well  known,  as  demonstrated  experimentally,  for  example 
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Table  2.  Adsorption  energies  E a(js  (kcal/mol)  for  NBs  and  WC  base 
pairs  adsorbed  on  SWCNTs  (B97-D/6-31G*)-  (Note:  values 

were  systematically  decreased  by  about  15-17  kcal  moP1  for  single 
point  calculations  (B97-D/Def2-TZVP).) 


NB 

C(6,  5) 

C(9,  1) 

C(8,  3) 

G 

-20.72 

-20.36 

-20.96 

A 

-16.99 

-16.79 

-16.66 

C 

-14.64 

-14.20 

-15.72 

T 

-15.18 

-14.59 

-14.85 

WC  NB  base  pair 

C(6,  5) 

C(9,  1) 

C(8,  3) 

AT 

-28.51 

-27.57 

-29.03 

GC 

-29.66 

-29.01 

-29.05 

for  adenine  adsorption  on  graphite  [25],  and  by  empirical 
force-field  simulations  ([26]  and  references  therein),  where 
the  potentials  were  extensively  parameterized  for  London 
dispersion  interactions.  We  note  that  for  the  SWCNT 
chiralities  considered,  the  distances  of  the  NBs’  carbon, 
nitrogen  and  oxygen  atoms  from  the  surface  were  3. 0-3. 2  A, 
3. 2-3. 3  A,  3. 3-3. 6  A,  respectively,  with  an  average  nearest 
neighbor  distance  of  ca.  3.2  A,  consistent  with  B97-D  results 
for  NBs  adsorbed  on  graphene  and  for  n -stacked  NBs  [27]. 
The  values  are  somewhat  larger  than  those  calculated  with 
the  M0-2X  [28]  functional  for  NBs  adsorbed  on  armchair 
SWCNTs  [9]  and  zigzag  C(10,  0)  [7],  presumably  because 
of  the  lack  of  an  explicit  dispersion  correction.  Notably, 
calculations  carried  out  with  the  PBE  functional  resulted 
in  non-stacked  NB  structures  on  the  SWCNT.  A  vertical 
positioning  for  the  NB  relative  to  the  SWCNT  surface  was 
obtained  for  T  on  C(6,  5),  a  tilted  positioning  for  C  on  C(6, 
5),  and  similarly  for  the  G  and  A  NBs  on  C(8,  3)  (figure  SI 
in  supplementary  data  available  at  stacks.iop.org/Nano/23/ 
165703/mmedia).  Interestingly,  a  different  behavior  for  the 
pyrimidine-based  NBs  (C  and  T)  and  purine-based  (G  and 
A)  structures  for  varying  chiralities  was  shown,  which  may 
influence  the  affinity  towards  self-stacking  in  water  [29];  this 
is,  however,  beyond  the  scope  of  this  study. 

Adsorption  energies  were  calculated  from  E-^  = 
£(SWCNT  -  NB)  -  //(SWCNT)  -  //(NB),  as  summarized 
in  table  2.  Results  for  C(5,  0),  for  comparison  with  previous 
work  but  here  using  B97-D,  were  calculated  as  —19.2,  —15.0, 
—  14.8  and  —13.8  kcalmoP1  for  G,  A,  T,  C,  respectively. 
The  trend  in  the  calculated  adsorption  energies  of  the  NBs 
on  same  diameter  C(6,  5)  and  C(9,  1),  and  on  C(5,  0),  was 
predicted  as  G  >  A  >  T  >  C,  consistent  with  previously 
calculated  values  for  graphene  and  C(5,  0)  [8,  9,  27],  The  //acjs 
values  were  lower  for  a  higher  SWCNT  curvature  because  of 
a  smaller  contribution  of  dispersion.  The  predicted  adsorption 
trend  is  consistent  with  experiments  on  the  interaction  of 
ssDNA  on  a  graphitic  surface  [30,  31].  Similarly,  the  C(8, 
3)  results  are  consistent  with  isothermal  titration  calorimetry 
values,  using  samples  of  graphene,  where  in  the  order 
of  G>A>C>T  the  positions  of  C  and  T  seemed 
to  be  interchangeable  [32].  Note  that  the  values  of  E^ 
were  systematically  decreased  by  about  15-17  kcal  mol-1 
for  single  point  B97-D/Def2-TZVP  calculations,  but  the 
adsorption  trend  did  not  change.  Furthermore,  our  qualitative 
agreement  for  C(6,  5)  with  the  results  of  Antony  et  al  [27]  for 


Table  3.  Adsorption  energies  £a( js  (kcal  mol  *)  of  NBs  adsorbed  on 
C(5,  0)  using  EFP2. 


G 

A 

T 

C 

Electrostatic 

-5.89 

-1.82 

-3.98 

-4.02 

Repulsion 

6.72 

4.09 

8.29 

5.68 

Polarization 

-3.17 

-0.64 

-3.20 

-3.49 

Dispersion 

-24.15 

-19.04 

-27.20 

-19.67 

which  larger  graphene  sheet  models  showed  the  same  trend, 
implies  that  this  relation  will  not  be  affected  by  an  increase 
of  the  already  relatively  large  chiral  SWCNT  model  systems. 
The  order  of  adsorption  strength  for  the  C  and  T  is  reversed  for 
C(8,  3)  compared  to  C(6,  5)  and  C(9,  1)  (table  2),  however  of 
a  slightly  different  diameter.  The  dispersion  contributions  for 
C(6,  5)  and  C(8,  3)  showed  a  similar  trend,  e.g.  for  C(6,  5): 
-27.50,  -25.10,  -22.70,  -20.50;  C(8,  3):  -28.20,  -23.00, 
—22.20,  —21.9  kcal  mol-1,  for  G,  A,  T,  C,  respectively.  The 
trend  can  be  explained  by  the  relative  polarizabilities  of  the 
NBs  for  a  given  carbon-based  nanostructure.  The  results  of 
the  relative  trend  are  indeed  consistent  with  the  trend  in  the 
experimentally  determined  ionization  potentials  of  the  NBs 
that  were  compiled  by  Roca-Sanjuan  et  al  [33],  This  trend 
may  also  be  related  to  the  so-called  7T -misalignment  [34]  in 
SWCNTs  of  varying  chirality. 

At  the  same  time,  recent  experimental  data  for  single¬ 
strand  homo-oligonucleotides,  d( A)  12,  <7(T)i2,  d(C)n  and 
d{ G)i2,  exhibited  a  different  trend  in  adsorption,  where  a 
stability  order  ofG>C>A>T  was  observed  [35].  This 
was  explained  in  terms  of  electrostatic  interactions  between 
the  phosphate  backbone  and  the  aqueous  environment,  or 
base  stacking  within  a  single  strand  due  to  changes  in  the 
SWCNTs  as  compared  to  a  flat  graphitic  surface;  this  is, 
however,  still  not  fully  elucidated.  It  is  interesting  to  point  out 
that  the  EFP2  results  (see  table  3)  demonstrated  a  similar  trend 
ofG>C>A>Tin  adsorption  energies  upon  exclusion 
of  dispersion  interactions  (which  may  not  be  accurately 
calculated  in  EFP  from  TDHF  [36]).  The  complexity  of  the 
theoretical  approximations  is  thus  highlighted,  and  we  also 
note  that,  although  adsorption  is  primarily  driven  by  tt-tt 
stacking,  the  experimental  environment  may  affect  the  order 
of  adsorption. 

The  enrichment  of  C(6,  5)  using  genomic  SaDNA  [5] 
demonstrated  an  important  advantage  practically,  and  there¬ 
fore  an  understanding  of  the  intrinsic  adsorption  energies 
could  provide  a  first  step  towards  further  exploration  of 
such  hybrid  materials.  Optimized  AT  and  GT  WC  base  pairs 
and  the  corresponding  structures  adsorbed  on  SWCNTs  are 
shown  in  figures  3(a)  and  (b),  respectively.  The  calculated 
and  experimental  geometries  [37]  of  the  WC  base  pairs 
are  consistent  with  our  results  (summarized  in  table  4). 
Upon  adsorption  on  the  SWCNTs,  curving  of  the  WC  base 
pairs  from  planarity  to  maximize  dispersion  interactions  is 
demonstrated  in  figure  3(b).  The  effects  of  adsorption  on  the 
hydrogen  bonding  in  the  WC  base  pair  are  small,  although 
slightly  larger  for  C(8,  3)  than  for  C(6,  5).  Note  that  the 
adsorption  energies  in  this  case  are  non-additive,  namely  the 
adsorption  energies  of  G  and  C  do  not  add  up  to  those  of 
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Table  4.  Intermolecular  distances  (A)  in  AT  and  GC  base  pairs  using  B97-D/6-31G*. 


AT  base  pair 


N6-04 

N1-N3 

H6-04 

N1-H3 

H2-02 

AT 

2.92 

2.82 

1.90 

1.77 

2.76 

Adsorbed  on  C(6,  5) 

2.92 

2.83 

1.91 

1.81 

2.68 

Adsorbed  on  C(8,  3) 

2.88 

2.75 

1.95 

1.73 

2.66 

GC  base  pair 

Ol-Nl 

N2-N2' 

N3-03 

Ol-Hl 

H2-N2' 

GC 

2.79 

2.92 

2.92 

1.74 

1.88 

Adsorbed  on  C(6,  5) 

2.81 

2.92 

2.88 

1.78 

1.91 

Adsorbed  on  C(8,  3) 

2.82 

2.91 

2.90 

1.77 

1.88 

AT  adsorbed  on  C(6,5) 


GC  adsorbed  onC(6,5) 


AT  adsorbed  on  C(8,3) 


GC  adsorbed  onC(8,3) 


Figure  3.  Hydrogen  bonding  for  (a)  WC  base  pairs  AT  and  GC,  and 
(b)  base  pairs  adsorbed  on  C(6,  5),  C(9,  1),  and  C(8,  3)  SWCNTs. 


a  GC  WC  base  pair.  This  is  due  to  the  hydrogen  bonding 
within  the  WC  base  pair  that  changes  the  electronic  structure 
and  strength  of  adsorption.  For  example,  the  Mulliken  partial 


atomic  charges  of  G  and  C  separately  adsorbing  on  C(8,  3) 
are  O.Ole,  while  the  respective  value  for  a  GC  WC  base  pair 
is  0.04e.  For  further  validation,  in  mimicking  experimental 
characteristics,  calculations  of  a  double  stack  were  also 
carried  out,  assumed  to  maintain  the  twist  of  the  double 
helix  in  the  dsDNA  (see  figure  S2  available  at  stacks.iop.org/ 
Nano/23/1 65703/mmedia).  The  double  base  pair  adsorption 
interaction  energies  demonstrated  the  same  trend  as  the 
single  WC  base  pairs  (representative  adsorption  energy  results 
are  shown  in  figure  S3  available  at  stacks.iop.org/Nano/23/ 
165703/mmedia),  indicating  that  the  adsorption  is  not  affected 
by  base  pairs  that  are  positioned  further  away  from  the 
SWCNT  surface.  However,  upon  possible  partial  unfolding 
of  the  double  helical  structure  in  a  realistic  environment 
even  stronger  adsorption  is  expected,  thus  motivating  future 
exploration  of  SWCNT  dispersion  with  genomic  salmon 
DNA. 

4.  Conclusion 

In  this  work  we  reported  calculations  of  NB  adsorption  on 
chiral  C(6,  5),  C(9,  1),  and  C(8,  3)  SWCNTs  to  elucidate 
trends  in  the  intrinsic  nonbonded  interactions.  Although 
the  geometries  of  the  hybrid  material  model  structures 
were  qualitatively  consistent  with  previous  work,  the 
importance  of  correcting  the  exchange-correlation  functional 
for  dispersion  was  clearly  demonstrated.  Limitations  in 
modeling  the  dispersion  interactions  by  considering  the 
SWCNT  as  a  molecular  model  may  mask  subtle  effects  in  a 
molecular-macroscopic  system. 

The  trend  in  calculated  adsorption  energies  of  the  NBs 
on  same  diameter  C(6,  5)  and  C(9,  1)  SWCNT  surfaces  was 
predicted  as  G  >  A  >  T  >  C,  which  is  consistent  with  related 
computations  and  experimental  work  on  graphitic  surfaces. 
The  calculated  dispersion  contributions  followed  this  trend, 
explained  by  the  relative  polarizability  of  the  NBs  for  a 
given  carbon-based  nanostructure.  The  order  of  adsorption 
strength  for  the  C  and  T  NBs  was  reversed  for  C(8,  3), 
which  is  difficult  to  assess  in  a  direct  comparison  because  of 
change  in  diameter.  However,  recent  experimental  data  on  the 
adsorption  of  single-strand  short  homo-oligonucleotides  on 
SWCNTs  demonstrated  a  trend  ofG>C>A>T  [35].  The 
possible  importance  of  electrostatic  interactions  was  partially 
captured  by  the  EFP2  results  when  excluding  the  dispersion 
contribution,  which  was  calculated  at  the  TDHF  level,  to  be 
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improved  upon  in  future  work.  Finally,  because  the  calculated 
adsorption  energies  for  Watson-Crick  base  pairs  have  shown 
little  effect  upon  adsorption  of  the  base  pair  farther  from 
the  surface,  the  results  on  SWCNT  sorting  by  SaDNA  could 
be  indicative  of  partial  unfolding  of  the  double  helix  upon 
adsorption  on  the  surface. 
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